NONLINEAR DIFFUSIVE-DISPERSIVE LIMITS FOR 
MULTIDIMENSIONAL CONSERVATION LAWS 



JOAQUIM M. CORREIA^ AND PHILIPPE G. LeFlOCH^'^ 

Abstract. We consider a class of multidimensional conservation laws with vanishing 
nonlinear diffusion and dispersion terms. Under a condition on the relative size of 
the diffusion and dispersion coefficients, we establish that the diffusive-dispersive 
solutions are uniformly bounded in a space {p arbitrary large, depending on the 
nonlinearity of the diffusion) and converge to the classical, entropy solution of the 
corresponding multidimensional, hyperbolic conservation law. Previous results were 
restricted to one-dimensional equations and specific spaces L^. Our proof is based 
on DiPerna's uniqueness theorem in the class of entropy measure- valued solutions. 



1. Introduction. 

Nonlinear hyperbolic conservation laws arise in the modeling of many problems 
from continuum mechanics, physics, chemistry, etc. The equations become par- 
abolic when additional dissipation mechanisms are taken into account: diffusion, 
heat conduction, capillarity in fluids, Hall effect in magnetohydrodynamics, etc. 
From a mathematical standpoint, hyperbolic equations admit discontinuous solu- 
tions while parabolic equations have smooth solutions. Discontinuous solutions, 
understood in the generalized sense of the distribution theory, are usually non- 
unique. It is therefore fundamental to understand which solutions are selected by a 
specific zero diffusion-dispersion limit. In this paper we address this issue for mul- 
tidimensional, scalar conservation laws, and review previous work on the subject 
restricted to one-dimensional equations. 

Consider the Cauchy problem 

dtU + dWf{u) ^0, (2;,t) e R'^ X IR+, (1.1a) 
u{x,0) ^uo{x), xem.'^, (1.1b) 

where the unknown function u is scalar- valued. Smooth solutions to (1.1) also 
satisfy an infinite list of additional conservation laws : 

dtrjiu) + div q{u) ^0, V„g = V„77VJ-, (1.2) 
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where 7] is a convex function of u. For discontinuous solutions, Kruzkov [5] shows 
that (1.2) should be replaced by the set of inequalities 

dtviu) + div q{u) <0, (1.3) 

which select physically meaningful, discontinuous solutions. The condition (1.3) is 
called an entropy inequality; it is motivated by the second law of thermodynamics, 
in the context of gas dynamics. By definition, an entropy solution of problem (1.1) 
satisfies (1.1) in the sense of distributions, and additionally (1.3) for any entropy 
pair (77, q) with convex t]. 

Consider the following approximation of (1.1) obtained by adding a nonlinear 
diffusion, b : IR"^ IR"^, and a linear dispersion to the right hand side of (1.1a) : 

dtu + divf{u)=diy(ebj{Vu) + 5dlu) , (x,t) e R'' x IR+, (1.4a) 

u{x, 0) = ul'\x), X G M'^. (1.4b) 

Let u^-^ : Wi'^ X [0, T] ^ IR be smooth solutions defined on an interval [0, T] with 
a uniform T independent of s,6. In (1.4b), Uq is an approximation of the initial 
condition Uq in (1.1b). 

Our main purpose is to derive conditions under which, as e > and S tend to 
zero, the solutions u'^'^ converge in a strong topology to the entropy solution of (1.1). 
When e = 0, equation (1.4a) is a generalized version of the well-known Korteweg- 
de Vries (KdV) equation, and the solutions become more and more oscillatory as 
(5^0: the approximate solutions do not converge to zero; see Lax and Levermore 
[6]. When 6 = 0, (1.4a) reduces to a nonlinear parabolic equation, resembling the 
pseudo- viscosity approximation of von Neumann and Richtmyer [8] : in that regime, 
the solution converges strongly to the entropy solution. Therefore, to ensure the 
convergence of the zero diffusion-dispersion approximation (1.4), it is necessary 
that diffusion dominate dispersion. Indeed the main result of the present paper 
establishes, under rather broad assumptions (see Section 3, Theorem 3.1-3.3), that 
the solution of (1.4) tends to the classical entropy solution of (1.1) when e,5 ^ 
with \S\ « £. 

For clarity, the main assumptions made in this paper are collected here. First 
concerning the fiux function we shall assume 

(ill) for some Ci, > and m > 0, \f{u)\ < Ci + C( iMp"^ for all m e H. 
For the diffusion term, we fix r > and assume 

{H2) for some C2 > 0, C2 |A1'^+^ < A • h{\) < C3 |Ar+^ for all A e W^. 
In the case < r < 2, we will need also 

(ifa) Dh{\) is a positive definite matrix uniformly in A G H''. 

We remark that the diffusion hj{Vu) = dxjU satisfies {Hz). 

The case d = 1 of one-dimensional c;quations was treated in the important pa- 
per by Schonbek [9], where, in particular, the concept of Young measures is 
introduced together with an extension of the compensated compactness method 
for conservation laws. We follow here LeFloch and Natalini [7] who, also for one- 
dimensional equations, developed another approach based on DiPerna's uniqueness 
theorem for entropy measure- valued solutions [2] (see Section 2 for a review). Specif- 
ically one uses a generalization of DiPerna's result to functions, due to Szepessy 
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[10]. The present paper therefore relies on a method of proof that was successful 
first in proving convergence of finite difference schemes : Szepessy ([10] and the 
references therein by Szepessy and co-authors) and Coquel and LeFloch [1]. 

Recent work by Hayes and LeFloch (see [3,4]) treats the transitional case where 
both terms, the diffusion and the dispersion, are in balance. Convergence results 
in this regime cannot be obtained by the measure-valued solutions approach. 

2. Entropy Measure- Valued Solutions. 

We include here, as background for further reference, basic material on Young 
measures and entropy measure-valued (e.m.-v.) solutions. First of all we will 
need Schonbek's representation theorem for the Young nicasm'es associated with a 
sequence of uniformly bounded in L^. The corresponding setting in was first 
established by Tartar [11]. In the whole of this subsection, q e (1, oo) and T < oo 
are fixed. 

Lemma 2.1. (See [9].) Let {uk} he a hounded sequence in L°°((0, T); L«(IR'')). 
Then there exists a suhsequence still denoted hy {uk} and a weakly * measurahle 
mapping v : M'* x (0, T) Prob (H) taking its values in the space of non-negative 
m,easures with unit total mass (probability measures) such that, for all functions 
g S C(]R) satisfying 

g{u) = 0{\u\^) as |u| — »■ 00, for some s G [0,q), (2.1) 

the following limit representation holds 

lim // g{uk(x,t))(t){x,t) dxdt = jj {v(^^^t),9) 4>(x,t) dxdt (2.2) 

's-^oo J JlRd X (0,T) JJIR<*X(0,T) 

for all(j)e L^iM!^ x (0,T)) n L°°(]R'^ x (0,T)). 

Conversely, given v , there exists a sequence {uk} satisfying the same conditions 
as above such that (2.2) holds for any g satisfying (2.1). 

We use the notation {v{x,t) i9) /]r fl'(^) di'(x,t) ) which therefore describes weak* 
— \img{uk)- "Weak ★ measurable" means that the real-valued function {y(x,t)ig) 
is measurable with respect to (.x,t), for each continuous g satisfying (2.1). The 
measure-valued function i'[x,t) is called a Young measure associated with the se- 
quence {uk}- The following result reveals the connection between the structure of 
y and the strong convergence of the subsequence. 

Lemma 2.2. Suppose that v is a Young measure associated with a sequence {uk}, 
bounded in L°°((0, T); L«(]R'')). For u G L°°{{0,T); Li{W^)), the following state- 
ments are equivalent: 

(i) limfe^oo Wfe =u in L^{{0,T); L'^^^ (IR-"')); for all s < oo and p G [l,q); 

(ii) z^(x,t) = Su{x,t) for a.e. {x,t) eW^ x (0,T). 

In (ii) above, the notation ^„(x,t) is used for the Dirac mass defined by 

{Su{x,t),9) = 9{u{x, t)) for all g G C{JR) satisfying (2.1). 

Following DiPerna [2] and Szepessy [10], we define the e.m.-v. solutions to the first 
order Cauchy problem (1.1). 
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Definition 2.1. Assume that / S C{JR)'^ satisfies the growth condition (2.1) arid 
uo € {JR'^) n L'l (SRf) . A Young measure v associated with a sequence {uk}, which 
is assumed to be bounded in L°°((0, T); L«(]R'')), is called an entropy measure-valued 
(e.m.-v.) solution to the problem (1.1) if 

9t(i^(.), |w-fc|) +div(i^(.), sgn(?i- fc)(/(M) -/(fc))) < for all keWi, (2.5a) 
in the sense of distributions on K'' x (0, T) and 

lim - / / (i^(x,s)) \u — uo{x) \ ) dxds = 0, /or a^/ compact set K C H''. 

(2.5b) 

A function m G L°°((0, T); L1(]R'') nL«(]R'^)) is an entropy weak solution to (1.1) 
in the sense of Kruzkov [5] and Volpert [12] if and only if the Dirac measure is an 
e.m.-v. solution. In the case q — +00, existence and uniqueness of such solutions 
was shown in [5]. The following results on e.m.-v. solutions were established in 
[10]: Proposition 2.3 states that e.m.-v. solutions are actually Kruzkov's solutions. 
Proposition 2.4 states that the problem has a unique solution in L'^. 

Proposition 2.3. Assume that f satisfies (2.1) and uq G L^(IR'') n L«(IR''). 
Suppose that v is an e.m.-v. solution to (1.1). Then there exists a function 
u e L°°((0, T); Li(]R'^) n L«(IR'^)) such that 

^{.,t) = K{x,t) for a. e. {x,t) & W x {Q^T). (2.6) 

Proposition 2.4. Assume that f satisfies (2.1) and uq G L^(IR'^) nL«(IR'^). Then 
there exists a unique entropy solution 

u e L°°((0, T); L^{JR'^) n L«(]R'^)) 

to (1.1) which, moreover, satisfies 

l|w(*)llip(iR<i) < ll^io||LP(iR<i) for a.e. t e (0,T) and all p £ [l,q\. (2.7) 

The measure-valued mapping V{x,t) — 5u{x,t) is the unique e.m.-v. solution of the 
same problem. 

Combining Propositions 2.3 and 2.4 and Lemma 2.2, we obtain the main con- 
vergence tool : 

Corollary 2.5. Assume that f satisfies (2.1) anduo G L'^ {W'') D L'' {M,'^) for q > 1. 
Let {uk} be a sequence bounded in L°°{{0, T); i^(IR'')) and let u be a Young measure 
associated with this sequence. If v is an e.m.-v. solution to (1.1), then 

\im.Uk = u m L*((0,r);Lfg^(IR'')) for all s < 00 and all p e [l,q), 
where u G L°°((0,T); L^{JR'^) n L«(]R'^)) is the unique entropy solution to (1.1). 
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3. Convergence Results. 

Throughout it is assumed uq E L^(iR'^)riL'^{iR'^) and the initial data in (1.4b) are 
smooth functions with compact support and are uniformly bounded in L^(]R'') fl 
L9(IR'^) for some q > 2. While in previous works [9, 7], a single value of q was 
treated, we can here handle arbitrary large values of q. For simplicity in the pre- 
sentation, we will always consider exponents q of the form 

q = 2 + n{r- 1), 

where n > is any integer. Therefore, when the diffusion is superlinear, in the 

sense that {H2) holds with r > 1, then arbitrary large values of q are obtained. 
Restricting attention to the diffusion-dominant regime 6 = 0{e), we suppose that 
Uq'^ approaches the initial condition uq of (1.1b) in the sense that : 

lim uf = Uq in L^{B.'^) f) L«(]R'^), 

^^°+ (3.1) 

Iko IIl2(IR.<') ^ IKolliSfipjd). 

The following convergence theorems concern a sequence u^'^ of smooth solutions to 
problem (1.4), defined on M'* x [0, T] and decaying rapidly at infinity. 

First consider the hypothesis {H2) with r > 2, that is the case of difi^usions with 
(at least) quadratic growth. 

Theorem 3.1. Suppose that the flux / satisfies {Hi) with m < q (which is al- 
ways possible when r > 1 by taking q large enough). Suppose that the diffusion 

b satisfies {H2) with r > 2. If 5 = o{e^), then the sequence u^'^ converges in 
((0,T); LP(5R'^)) , for all s < 00 and p < q, to a function 

uEL°° ((0, T); L\Wi'^) n L«(IR'^)) , 

which is the unique entropy solution to (1.1). 

Observe that m and q can be arbitary large in Theorem 3.1. To treat the case 
r < 2, we need the additional condition (H3) on the diffusion. First for difl[Usion 
with linear growth (r = 1). we obtain a result in the space LF' : 

Theorem 3.2. Suppose that f satisfies {Hi) with m < 1, and b satisfies {H2)-{H^) 
with r = 1. If 6 = o(£^), then the sequence u^'^ converges in ((0, T); L^'(]R'*)), 
for all s < 00 and p <2, to a function 

u€L^ ((0, T):L\TR^) n L''{Wi'^)) , 

which is the unique entropy solution to (1.1). 

In particular Theorem 3.2 covers the interesting case of a linear diffusion and 
a linear dispersion with an (at most) linear flux at infinity. More generally, for 
general r > 1 we establish that : 

Theorem 3.3. Suppose that f satisfies {Hi) with m < < q, and b satisfies 

{H2)-{Hz) for some r > 1. If 5 ~ o(e~) then the sequence u^'^ converges in 
L" ((0, T); LP(]R'^)), for all s < 00 and p < q, to a function 

u€L^ ((0, T):L\-iR'^) n L«(1R'^)) , 
which is the unique entropy solution to (1.1). 

Our method of proof can also be extended to a general diffusion b = b{u, Vu, D^u). 
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4. Convergence Proofs. 

The superscripts e and 5 are omitted in this section, except when necessary. In 
the proof, we make frequent use of the fohowing computation. Multiply (1.4a) by 
r]'{u) where : IR — > IR is a sufficiently smooth function and define g : IR ^ IR'* by 
q'j = 77' f'y We have 

dtr){u) = - 77'(u) div/(u) + e^a,/r/'(u) 6,(Vw)) - d,Ji'{u) b,{\Iu) 

j 

= - div q{u) + £ ^ d^lr)\u) bj{S/u)) - e r)"{u) ^ d^.u bj{Vu) 
j 



\2 



thus 



dtr]{u) + divg(u) =e div(?7'(w) 6(Vu)) — er]"{u) Wu ■ 6(Vu) 

- ^E^"W d.,{d.,uf + Sj2d.,{v'{u) din) , ('^•1^) 
J j 

The last two terms in the right hand side of (4.1a) take also the form 

^ 5^ V"'{u) {d.,uf - 3 d.,{v"{u) {d.,u) + 2 d^v\u) d.,u) . (4.1b) 



When 1] is convex, the term containing r]"{u) has a favorable sign : the diffusion 
dissipates the entropy rj. 

We begin by collecting fundamental energy estimates in several lemma. 



Lemma 4.1. Let a > 1 be any real. Any solution of (1.4a) satisfies, for t G [0, T], 
I ^-^^^X^dx + ae f f \u\°'-^Wu- b{Wu)dxds 



For a>2, the last term in the above identity can be replaced by 
a {a -I) 



sf f sgn(u) V (a^^w)^ rfa;rfs. (4.2b) 

Jo JtR'' ^ 
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Proof. Integrate (4.1a) over the whole of IR"^ with r]{u) = ■'^qij- : 
at J]Rd a + 1 

= -ae f \u\"-^Wu ■ b{Wu) dx - ^ f V \u\''-^d^p^^u)^ dx, 

which yields (4.2a) after integration over [0,t]. One may use (4.1b), instead, to 
obtain (4.2b). □ 

Choosing a = 1 in Lemma 4.1, we deduce immediately a uniform bound for u in 
L°°((0, T);L^{W^)) together with a control for both Vm ■ b{Vu) in L\W^ x (0, T)) 
and Vu in L'-+\TR'^ x (0,T)). 

Proposition 4.2. For any solution of (1.4a) and t G [0,T], we have 

I u{tfdx + 2e[ I Vu ■ b{Vu) dxds = [ u^dx (4.3) 
and, assuming {H2), 



\Vu\''+^dxds < C 



JTR'i 



dx. 



(4.4) 



To derive additional a priori estimates, we use another value of a, motivated by 
controling the dispersive term in (4.2b) with Holder inequality, as follows : 



nsgn(u) 2 Y] i^xju)^ dxds 

< / / 1"!""^ |Vt<|^ dxds 
Jo Jtr'' 

Jo Jm.'' i Uo Jtb 



(4.5) 



Wul^P' dxds 



To take advantage of (4.4), we can choose 3p' = r + 1 provided r > 2. Then 
p = ^3^, so {a — 2)p = (r + 1)^5|- Therefore it is rather natural to take the 
exponent a = r for the entropy, where r is given by the diffusion term. Thus we 
deduce from Lemma 4.1 a natural estimate for |u(f)|''+^, involving the combination 
5 of 6 and e. 

Proposition 4.3. Assume that {H2) holds with r > 2 and uq € L^~^^{5R'^). For 
t G [0,T], we have 

I |u(t)r+^dx + (r + l)r£ / / \uy-^Vu-h{Vu)dxds 
Jtr'' Jo Jtr'' 

< Ci(wo) (^l + <5e"^ maxjl, t Ci(mo) (l + (5 e^^^ll ~ 1^ (4-6) 
:= Hi(^S £-4t) 
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and ^ 

ell \u\'-^\Vu\''+^ dxds<—^—Hi(6e-^) , (4.7) 
Jo JjR" (r + 1) r V / 

where C > is some fixed constant and 

In particular Proposition 4.3 shows that, if uq G L"^ H 1/'"+^ and d = O(e^), 
then u{t) e L'^+i uniformly for alH > 0. 

To motivate the forthcoming derivation, let us consider the special case r = 2. 
Then (4.6) gives us an estimate. Returning to the original inequality (4.5), but 
now with the new value a = 3, we now can estimate the dispersive term in (4.2b) 
directly in view of the estimate (4.7). In this fashion, we deduce an estimate 
from Lemma 4.1. This argument can be continued inductively to reach any space 
L". 

Actually Propositions 4.2 and 4.3 are the first two cases of a general result derived 
now. We define, for n > 1, 

Ho(d £"^) = Co(uo) := ||wo||i2(]Rd), 

Hn[d s-^y.= C„{uo) (l + 6 max 1^1, t C„(uo) (l + (5 e"^) ^| 

and 

r, r ^ Jii i|n(r-i)+2 n(r - 1) + 2 n(r -1) + 1 
Cn[uo) :=max<^ lluo|li„(.-i)+2(]Rd) 



[(n - l)(r - 1) + 2]^ [{n - l)(r - 1) + 1]^ 

(4.8) 

Here C > is some fixed constant. Note that iJ„ and C„ are uniformly bounded 
in e, S provided uq e L'^ Ci L'^(r-i)+2 ^^^^ ^ ^ O(e^). 

Proposition 4.4. Assume that {H^) holds with r > 2 and uq G L'^{M.'^). For 
t G [0, T] and n > such that n{r — 1) + 2 < we have 

+ e(n(r-l) + 2)(n(r-l) + l) / / \u\''<^'-'^^ ■ b{\7u) dxds (4-9) 



<Hn{5 



Jtr'' 



and 



iwr^'-i) ivur+' dxds 

/o Jm^ (4.10) 
< C (n(r - 1) + 2)-\n{r - 1) + l)-^if„(<5 e"^) . 
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Proof of Propositions 4.3 and 4.4. Note first that (4.10) is an immediate 
consequence of (4.9) and the hypothesis (i?2)- If n = 0, (4.9) coincides with (4.3) 
in Proposition 4.2. For n = 1, the estimate is Proposition 4.3. 
To estimate the term in (4.2b), with a = r, we use (4.5) : 



/ |u(t)r+^dx + (r+l)re / / \u\''-'^ Vu ■ b{Vu) dxds 

JlR<* Jo Jm'* 



< [ iwor+^rfx 

^ (r + l)r(r-l) ^ 



(4.11) 



JTR'' 



Wu\'^~^^ dxds 



JTR'' 



3 
r+1 



By {H2) the second term in the left hand side of (4.11) is positive. Integrate (4.11) 

over [0, t] and use (4.4) : 

ll"lll,^+i(IR<ix(0,T)) 

<t\\uo\\ltl.^^.^ 



jr + l)r-(r - 1) , 2 

2 * ['^ ll"0|lL2(IR'i) 



t (C|l«0|li2(iRd))'"'' SS rll ||u||2.^1(IRdx(0,T)) 



+ 

< t C,{uo) (l + 5 e-^ (lKII2t+i(H^x(o,T))) 
Observe that, for X > 0, the inequaUty 

X < K (1 + AX"^^ . 
where < 9 < r + 1 and K > 0, implies 

X < maxjl, [ii'(l + A)]^^T:^|. 



(4.12) 



Thus we deduce 



I"IIl^+i(IR<*x(0,T)) - 



|l, [tCi{uo)(l + Se-^)] ' I 



and, returning to (4.11): 

/ \u{t)\''+Ux + {r + l)re 
< Ci(uo) ( 1 + max < 1 



/ / \u\''-^Vu-b{\7u)dxds 
Jo Jtr-' 

tCx{uo) (l + ^e-^i) 



This completes the proof of (4.6). 
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This argument can be iterated. We return to the dispersive term and make an 
estimate similar to (4.5), but now having in view to apply (4.10), already established 
for n = 1 : 



nsgn(u)|u|" ^ y2 {dx,uf dxds 

f f |Vu|^ dxds 

Jo Jtr'' 



< 



< 



(4.13) 



77 

Jo Jm'^ 



dxds 



[ [ \u\^P' \Vu\^P' dxds 
Jo JtR'' 



where we choose 3p' = r + 1 and 7p' = r — 1, so {a — 2 — ^)p = — 2 — 3 ^rpx) r^i • 
Then (4.2b) gives 

/ \u{t)\''+^dx + {a + l)ae [ [ \u\"-^ Wu ■ b{Wu) dxds 

ilR<i Jo JtR'' 

< [ \uo\"+'^dx 

1 — 2 

(a + l)a(a - 1) / / r. a 



+ 



2 [(r + l)r] '■+1 



"'IR'' 



|(a-2-7)p^2,^^ 



We choose a so that a + l = (a — 2 — 'y)p , i.e., a = 2j — 1 . 
Integrating (4.14) over the interval [0,t], we obtain 



(4.14) 



ll'"llL2r(]Rdx(0,T)) 



+ 



L2'-(IRd) 

r(2r- l)(2r-2) 
[(r + l)r]^ 



< tC2{uo)(l + 6e 4i (||M||^^2.(K,^(o_r))) 



with ^2(^*0) max^ \\uo\\%.,^.., '^''"^^ ^^"^^ (ci/i((5 



_3_ \ \ r+1 

r+1 



By (4.12), we obtain again 



II l|2^ 

ll^llL2r(]Rcix(0,T)) 

Then (4.14) gives 



< max < 1 , 



tC2{uo) (l + ^e-^) 



r+1 
3 



f \u{t)f''dx + 2r{2r-l)e f [ \uf^''-'^^ Wu ■ b{Vu) dxds 
Jtr'' Jo Jtr'' 

t C2{uo) (1 + 5 



< C2(uo) 1^1 + (5 £ '■+1 max <{ 1, 
:= H2(5 £-4t) . 



e '■+1 
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This establishes (4.9) for n = 2. The general case follows by induction on n. □ 

We are now concerned with the case where the diffusion exponent in {H2) satisfies 
r < 2. In this situation, we require the assumption (Hs), which for instance is 

satisfied by &j(V'u) = d^ -u. 

Proposition 4.5. Suppose that {Hi)-[H3) hold with m and r such that m < 
and r > 1. For t e [0, T], we have 



2 , 2(r+2) 

dx +£ -+1 



T 

Jm.<^ 



D^u dxdt < C, 



(4.15) 



+ 1 r+3 



(4.16) 



eTTT / |Vu(i)| 

/ \u{t)\^+'^ dx + e [ I \u\'^\Wu\''+^dxdt<c(l+ 5^-^ e 

Proof. We differentiate (1.4a) with respect to the space variable x : 

dtWu + div {f\u).Wu) = £V ^ (6j(V«)) + 8 (V«) 

3 j 

and, then, we multiply it by Vu and integrate over IR*^. After further integration 
by parts, we obtain 

1 d 



2dt 



■/ \\/u{t)fdx-[ Am/'(w) • Vuda; 

Thus, integrating on [0,t] using (iJi) yields 

/ \SIu{t)f dx + 2e / YVda:^u ■ Db{Vu) ■ Vd^^^udx 

_ f j \D''u\ lup-' |Vu| dxdt 
Jo Jtr'' 



< [ IVuor dx + 2Ci 



< 



|Vuo|^ dx 



and so, using (iJs), 



£ Jo JjR'' 



|u|^™"^|Vwl^ dxdt + 



\D u\ dxdt, 



JIR<* 



/ \Vuit)f dx +C5£ [ I \D^u\^ dxdt 
Jtr'' Jo Jtr'' 

< [ IVuol' dx + -f[ 
Jtsj^ ^ Jo Jtev^ 



lur™-^ ivur dxdt. 



By Holder inequality and for m < ^rpj 



/ |Vu(OI^ dx+C^e j j \D' 
Jtr'^ Jo Jm'' 



u\ dxdt 



< [ \Vuof dx + Ce-i 



f [ |Vwr+' dxdt] Iff \u 
Jo Jm'' i Uo JtR'^ 



dxdt 
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and now (4.15) follows from (4.3)-(4.4). 

To establish (4.16)-(4.17) we use (4.2a) for a > 1 : 



TR'' 



u{t)\''+' dx + C s / \u\''-'\Vu)\''+'dxdt 



JTR.'' 



<[ \uo\"+'^ dx + C 6 f f \u\"-'^Vu\\D\\dxdt. 
Jtr'' Jo i]R<* 

We evaluate the last term using {H2): 
\ur-^Vu\ \D^u\dxdt 



< 5 



10 Jm'^ 

rt 



{r + l)S 



r + 1 V C2 £ 
t 



< 



ill \u\"-^ \Vu\'^+^ dxdt +C"6'^£--r f f lul^-^lD^ul"^ 
2 Jo Jm'^ Jo Jtr'' 



dxdt 
dxdt. 



So, we have 



[ \u{t)\"+'^ dx + C e I I 
Jtr'^ Jo Jtr"^ 



u\"-^ |VMr+' dxdt 



< f \uo\''+Ux + CS'^e-^ [ f \u\"-^\D^u\^ 
Jtr'' Jo Jm'^ 

Taking a = 1 + we deduce 

/ \u{t)\'^+'^dx + Ce [ [ lu]"^ \Vu\''+'^ dxdt 
Jtr'' Jo Jtr'' 



dxdt. 





2r 




I I \u\'^dxdt 




I I \D'^u^ dxdt 


Jo Jtr"^ 




Jo Jm'' 



< I ^0^+^^ dx + C5re- 



The conclusion follows now easily. □ 

Proof of Theorems 3.1. We first derive (2.5a), based on the conservation law 
(4.1b) with a arbitrary convex function, ry, where we assume r]',r]",r]"' bounded 
functions on IR. We claim that there exists a bounded measure /x < such that 

dtv{u) + div q{u) — > IJ, in I?'(]R'' x (0,T)). 

Prom (4.1b), we obtain 

dtri{u) + divg(w) = e div{rj'{u) 6(Vu)) — e'i]"{u) Vu • 6(Vu) 
6 
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with obvious notation. For each positive 9 G C^iM!^ x (0,r)) we evaluate {fii,0) 
for i = 1,2,3. To treat /xi, we use Holder inequality with the exponent In 
view of {H2) and (4.4) of Proposition 4.2 and assumption (3.1), we get 

Kmi, ^)| < e / / V \v'iu) 6,(Vu) a^/l dxdt 
Jo Jtr-' j 

' f I |V6i| |6(Vu)| dxdt 
Jo Jm'' 



< Cs 



so 



i(IR<*x(0,T)) 

< C e~ ||V0||Lr+i(]Rdx(o,T)) 
For fi2, we use {H2) and the convexity of r] : 



|Vur+^ dxdt 



supp 



(^2, 6*) = -£ / / V 61 r]"{u) Vu • 6(Vw) dxdt < 0. 

Jo JlR'i 



For ^3, we use again Holder inequality, as follows 

\{f^3,0)\ 



< 



< 



\f I E k + ^n"{u) {d^,uf d^p + 2r?'(«) dip 

^ Jo Jm'' j I 

cs [ [ e |Vu|^ da;(ii+ CS [ [ JZldx^uf \d^p\ dxdt 
Jo Jtr'' Jo Jtr'' j 

Jo JTR'' A 



dxdt 



+ cs 



dxdt. 



so 



|(M3,^)| < CS\\6\\ r+l 



'// ivnr 

supp ^ 



dxdt 



+ cs\\^e\\ r+i 



supp 



+ CS 



r+l 



JlR<i 



dxdt 



[[ iv^r 

.J J supp^ 



dxdt 



Therefore, we conclude that 
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Finally the condition i5 = o{e~) is sufficient to imply the desired conclusion. 

Using a standard regularization of sgn(u) and \u — k\ (for fc G ]R), which fullfil 
the growth condition (2.1), we apply the limit representation (2.2) and conclude 
that ly satisfies (2.5a). 

To show (2.5b) we follow DiPerna [2] and Szepessy [10] 's arguments. We have 
to check that, for each compact K of Mf^, 

^^^7 I {'^{x.s), W ~ "^oix)]) dxds 

= hm lim - [ [ \u'^'^(x,s) — uq(x)\ dxds = 0. 
By Jensen's inequality, where m{K) stands for Lebesgue measure of K, we have 

- / / \u'^'^ [x, s) — uq{x)\ dxds 
t Jo Jk 



We will establish that 



lim lim — 

«^0+e^0+ t 



Jk 



{x, s) — Uq{x)) dxds j 



1/2 



Jk 



{x, s) — uo{x)) dxds = 0. 



Let Ki C ifj+i {i = 0,1, ... ) be an increasing sequence of compact sets such that 
Kq = K and Ui>oKi = H"^. We use the identity — u^ — 2uo{u — uq) = {u — uq)^ : 

- [ [ {u^'^{-,s) — uo)^ dxds 
t Jo Jk 

\u^'^{; s)\^ J uldx-2 J Uo {u^'^i-, s) - Uo) dx^ ds 



< 



\JK, 



< 



[ u^dx +- I [ Uo {u^'^{-, s) — Uo) dx 

JTR'^XKi t Jo JKi 



ds 



dx = 0, 



for all i = 0,1, where we used (4.3)-(3.1). 
Since 

lim / Uq 

Jm.'i\Ki 

we only consider the last term above. Take {dn}neTN C Cq°(]R'^) such that 

lim 9ri = Uo in L'^{JR'^), 

n— *oo 

Cauchy-Schwarz inequality gives 

/ Uo {u^'^{-, s) - Uo) dx 

JKi 



< 



'n\ |u^'''(-,s) -Uo\ dx 



Ki 



< \\U0 - 6'nl|L2(IR<i) (||w^''^(-, s)||L2(iRd) + ||wo||L2(]Rd)) 



e,5 



+ l|c'n||i2(iRd)||Wo' - Wo||L2(iRd) + 



fl 

Jo JKi 



On dsu'^'^ dxdr 
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In view of (4.3) and (3.1) 

||W0 - (||w^'''(-,s)||l,2(]R<i) + ||wo||i,2(]R^<i)) 

< 2||Uo||l,2(]R,ci)||UO — ^Tl||z,2(]R<i), 



which tends to zero when n ^ oo, and since hm£^o+ ll^o'' ~ '"o||i,2(]R<i) = by 
(3.1), it remains only to see that 



hm hm - 



Jo JKi 



d.,u^'^ dxdr 



0. 



We have, by (1.4a), 

9n dgU dxdr 



JKi 



JKi 



9n{ — div f{u) + s div b — S d^. u) dxdr 



I [ (ven-f{u)-eV0n-b + Sy"dl0n)udxdT 

Jo JKi V ' 



:= /ii +M2 + 

To deal with /ui, we use Holder inequality and (Hi) 



\V0n\ \f{u)\ dxdr 



JKi 



< C 



IVej dxdr + C 



Jk, 



Jk, 



iV^nl'-™ dxdr 



fl 

Jo JKi 



dxdr 



For fj.2, using {H2) and once more Holder inequality with(4.4)-(3.1), we get 
|V6'„| |6| dxdr 

\V0n\ IVwr dxdr 



'0 J 

<Cse 



11 

Jo JKi 



<C.e 



|V6'„I''+^ dxdT 



JKi 



.Jo JKi 



dxdr 



< Ce^ rll sril ||V6'„||ir+l(]R<i). 

Finally, for we use Cauchy-Schwarz inequality with (4.3)-(3.1) : 



'/7 

Jo JKi 



dxdr 



< 6 



\u\ dxdr 



JKi 



JKi 



dx dr 



<SS ||V^6„||l,2(]J^<i) ||Uo||l,2(]Rd), 
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thus 



lim - 

£^0+ t 



9r, dM^'^ dxdr 



JKi 



ds=:X 



+ c 
+ c 

and so 



q — m 



lL5^(IR<i)ll lli«(IR<*x(0,T)) 



r + 1 



(IR<ix(0,T)) 



where we have used (4.9) in Proposition 4.4. The desired conclusion when t —> 0+ 

follows. □ 

Proof of Theorems 3.2 and 3.3. In the previous proof, to establish (2.5a) we 
started with the identity (4.1) and the condition 6 = o{e^) as required, inmpar- 
ticular to control the term in (4.1b). We now keep the form (4.1a) instead (4.1b) : 
the terms t/ii and /i2 introduced in the previous proof do not change. We only 
need discuss /X3. It has now the form: 



The first term is bounded as follows 
/ Z_j ^ ^" (^) (^xjU u dxdt 



<C5 [ [ e\yu\ \D'^u\ dxdt 

Jo Jtr'' 

<CS [ [ n\D'^u\^ + -{9\yu\f dxdt 

Jo Jtr.'' ^^ 

( n + 2 1 2 

<Cb\\X£ ~ ^ £ --+1 



using (4.15) and (4.4), and we take [i = e and b = o{s'-+^ ). 
The second term in /U3 behaves better: 



9j2dxAv'iy')dlu) dxdt\ =:Y 



JTR'' 



< 6 



7 

Jm.'' 



< C5 



I I y,dl9ii'{u)d,^udxdt 

Jo JTR.'i ^ 

j j \Vu\\D^0\ dxdt + C 6 [ [ \Vu\^ \V9\ dxdt 
Jo Jtr,'' Jo Jtr'' 



'Y^dxflr]"{u) {dxjuf dxdt 
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thus 



r[ \yu\''+^ dxdt ^ +C5 Cj |Vur+' 

Jo JTRji Jo JtR'^ 



Y <CS 

<C5e~^. 

This completes the proof of Theorems 3.2 and 3.3. □ 



dxdt 
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